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ABSTRACT 

We derive the propagator for a massless, minimally coupled scalar on a .D- dimensional, 
spatially fiat, homogeneous and isotropic background with arbitrary constant deceleration 
parameter. Our construction uses the operator formalism, by integrating the Fourier mode 
sum. We give special attention to infrared corrections from the nonzero lower limit associated 
with working on finite spatial sections. These corrections eliminate infrared divergences 
that would otherwise be incorrectly treated by dimensional regularization, resulting in off- 
coincidence divergences for those special values of the deceleration parameter at which the 
infrared divergence is logarithmic. As an application we compute the expectation value of 
the scalar stress-energy tensor. 
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1 Introduction 



On the largest scales, the visible universe is well described by a homogeneous, isotropic and 
spatially flat geometry whose invariant element takes the form, 

ds 2 = g^dx^dx" = —dt 2 + a 2 {t)dx ■ dx . (1) 

Two derivatives of the scale factor a(t) have great importance, the Hubble parameter H(t) 
and the deceleration parameter q(t), 

H(t) = - , q (t) = -l-lL = -l + e{t). (2) 

Although the Hubble parameter may have changed by as much as 57 orders of magnitude 
from primordial inflation to now, the deceleration parameter is only thought to have varied 
from nearly —1, during primordial inflation, to nearly +1, during the epoch of radiation 
domination. And q(t) has been approximately constant for vast periods of cosmological 
evolution. It is therefore interesting to consider physics during an epoch of constant q(t). 

A field of great importance for primordial inflation is the massless, minimally coupled 
scalar, 

£ = -\d^d vV g^^-g = l<?a D - 1 - V" 3 . (3) 

This is true both as an approximation to the inflaton, on account of the very flat potential, 
and also because dynamical gravitons have the same kinetic operator [TJ. Note that we 
are assuming the spatial manifold is (D — 1) -dimensional in order to facilitate the use of 
dimensional regularization. 

The scalar propagator obeys the equation, 

^/^gDiA(x;x') = dJ y / = gg^dJA{x;^)\ = i5 D (x-x') . (4) 

An elegant solution to this equation has recently been obtained for the case of arbitrary 
constant e = —H/H 2 , which of course means constant deceleration [2J. This solution has 
been used to compute the one loop correction to the effective field equations in Einstein + 
Scalar [31 H]. Unfortunately, the solution diverges at certain values of e, even off coincidence 
and with the dimensional regularization in effect. The purpose of this paper is to elucidate 
the physical basis of this divergence and to provide a very simple correction for it. 

It is important to recognize that the propagator equation (j4j) only fixes iA(x;x') up to 
the addition of a homogeneous solution. There is no guarantee that a particular solution 
can be interpreted as the expectation value of the time-ordered product of <p(x) x <p(x') in 



1 



the presence of any state. One can see this even from the simple harmonic oscillator of 
nonrelativistic quantum mechanics [5]. The Heisenberg position operator can be expressed 
in terms of the initial position and momentum, 

Q(t) — Qo cos(ujt) + sin(u>t) . (5) 
mw 

The associated propagator in the presence of an arbitrary state can be expressed in terms 
of three real numbers, 

{i>\T[q(t)q(t')M) = sin(u\t-t'\) 

+a cos(ut) cos(ut') + j3 sin[u(t+t')] + 7 sin(a;t) sin(a;t') . (6) 

The parameters a, (3 and 7 are, 

<^={^\<&W) > P = o (^koPo+Po?o|^) , 7 = —5-0 (V'bol^) • (7) 

These parameters are not arbitrary. For example, the Uncertainty Principle implies, 

However, the right hand side of ([6]) solves the propagator equation for any choice of a, f3 
and 7, 

-m(J^ + u; 2 )iA (*;*') =i5(t-t') . (9) 

We seek a true propagator, rather than just a Green's function that solves (jl]). To ensure 
this we construct iA(x;x') from its operator mode sum in Section 2, assuming that the 
spatial manifold is M^" 1 . This gives precisely the result obtained previously In Section 
3 we note that the infinite space mode sum has an infrared divergence for any value of 
e < 2(D—1)/D. In most cases this is a power-law divergence that the automatic subtraction 
of dimensional regularization simply but erroneously discards. For certain special values of e 
the divergence becomes logarithmic, in which case the solution is infinite, even off coincidence 
and for general D. In Section 4 we solve the problem by formulating the mode sum on a finite 
spatial manifold. This changes the propagator by the addition of a series of homogeneous 
solutions which cancels the divergences. Section 5 uses the corrected propagator to compute 
the expectation value of the scalar stress tensor. Our discussion comprises Section 6. 
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2 Infinite Space Mode Sum 



The purpose of this section is to construct iA(x; x') from the canonical operator formalism, 
assuming the spatial manifold is We begin by working out the Hubble parameter 

and scale factor as functions of co-moving and conformal time. Then we give the free field 
expansion of <f(x). Taking the expectation value of the time-ordered product of (p(x) x <p(x') 
gives a Fourier integral expression for the propagator which we are able to evaluate by 
analytically continuing some integral identities. The section closes by noting the special 
values of e for which the resulting propagator diverges, away from coincidence and with 
dimensional regularization in effect. 

It is simple to reconstruct the Hubble parameter and the scale factor when e is constant. 
If we define H = H(0) and a(0) = 1 then relation (j2j) implies, 

IT 1 

H(t) = - °— and a(t) = l+eH t £ . (10) 

1 + eH t L 

The conformal time rj is defined by drj = dt/a(t). We can choose the zero of conformal time 
so that the following relations apply, 

H(v) = — H ° —. and a(r ] ) = — J , (11) 

For < e < 1 the universe is accelerating and conformal time approaches zero from below; 
for e > 1 the universe is decelerating and rj is positive with our convention. A very useful 
relation is, 

(l-e)Ha = -~. (12) 
V 

The massless, minimally coupled scalar on R D_1 can be expressed as a Fourier integral 
of plane waves, 

Voo ( V ,x) = J {u( V ,k)j n - s a(k) +u*( V ,k)e- in -*a\k)} . (13) 

Since (p^ is a minimally coupled massless scalar, it obeys the Klein Gordon equation d^oc = 
0. The mode functions are quite complicated for general, time-dependent e(t) [6]. However, 
for constant e they obey the simple equation 

D 1 , n _ D-1 



d\ -e + — ^ j (aT-y,, k) ) = 0, y= , i I I) 
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where k = \\k\\. The solution to this equation takes in this case a simple form 



■it 



^,k) = ^^a 1 -^Hi 1 \k\ V \). 
The creation and annihilation operators are canonically normalized, 

~a{k),a\k')] = (2tt) d - 1 5 D - 1 (k-k') . 



(15) 



(16) 



The state |0) which is annihilated by a; (A;) is known as Bunch-Davies vacuum |7j. It corre- 
sponds to the Heisenberg state of minimum excitation in the distant past. 
The infinite space propagator is easy to write as a mode sum, 



iA^ix^x') 



n 



9(t-t')ip 00 (x)ip 00 (x') + ^t'-tVooOc'VooOz) 
d D ~ X k 



n), 



D ik-(x—x') 



x ^(Ar ? )i/«(fc|r / |)^ 1 )(fc|V|)*+^(-A^)i7( 1 )(A;|7 7 |)*i7( 1 )(fc|r / '|) I, 



(17) 



where Arj = rj — rf . Now recall the D = 4 angular integral, 



d 3 k 
(2^p 



Sk- Ax 



/(Pll) 



2vr 2 



2 sm(fcAx) . . 
kAx K ' 



(19) 



Here and henceforth we define Ax = — Generalizing to D spacetime dimensions we 
have d D ~ l k = k D ~ 2 dk dQo-2, k = \\k\\, where 

dQ D ^ 2 = sm D - 3 (9 D _ 3 )d9 D - 3 sm D -\9 D ^)d9 D ^ ...d(f>, (20) 

where 9 D _ 3 , 9 D _ i , . . ., and are the angles on the sphere E> D ~ 2 . Making use of 

2tt^ 2(4 7 r)f- 1 r(f) 



dtt 



D-2 



T(^) T(D-V 



and Eq. (8.411.7) in [8] gives 



ik-Ax 



{2tt) d - 1 



f(\\k\\) 



2 D-2 



D-l 
71 2 Jo 



° Jd-3 (kAx) 

dkk D - 2 -f-—^f(k). 

(\kAx) — 



(21) 



(22) 



Hence the scalar propagator can be reduced to a single integral, 

D 



/;///Maa/)1 " ¥ r dkk D - 2 J2 ^ {kAx 



D-3 
2 



x |^(^)^T I S 1 )(A;|7 7 |)^ I S 1 >(A ; |VI)*+ ^(-^r 7 )^ 1 )(A ; |r 7 |)*^r j S 1 >(A;|VI) [■ 
To reach the final form we make the change of variable 



z = y r]7]' k - 
Recalling that 1 / (rja) = 

(rjrj 1 aa') 1 ^^ 



k/\x = , z 



y/r/rj 

l — e)H gives, 



k\r}\ 



k\v'\ 



2 D n^ Jo 



dzz°- 2 x 



D- 
2 
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Ax 






T „ 








D_ 




Ax 


*) 




2 A 


/vv' 





Arj) x I conjugate 



'l-e) 2 HH' 



--1 



(4tt; 



7ri2V f°° , n=i /Ax 

X rr^r X / dZ Z 2 X J c-3 

2 Vw 



D-3 

Ax \ 2 ^0 



x|^(Ar ? )i/( 1 )(y^^if( 1 )(^^* + ^(-Ar ? ) x (conjugate) 



(23) 



(24) 



(25) 



(26) 



By analytic continuation of relation (6.578.10) and employing relations (8.702), (8.407.1), 
(8.476.10) and (9.131.1) in [8] one can show, 

r(fji+i+u)r(fi+i-v) 



dx x" +1 J M (cx)#W (ax)H^ (bx)* 



r(/i+ 



3 (a+6) 2 -c 2 



(aby 



-2-P1 AH-l+f, /i+l-i/; 



4a6 



(27) 



Now make the following assignments for the various parameters in (!27j) . 



With these assignments we have, 
(a+b) 2 -c 2 1 



c 2 -{a-bf 



1 - 



Ax 2 - Arf 



(29) 



4a6 V Aab ) V 47/?/ 

Hence the mode sum ( l26l) for the propagator can be given the simple spacetime expression, 

Vi 



'l-e) 2 HH' 



r 



D-l 



Mr 



D-l 



(4tt) 



r(f; 



2^1 



D-l D-l D y 

Yv, v\ — ; 1 

2 '2 ' 2 ' 4 



Here and henceforth we define the quantity y = y(x;x') as, 



y{x;x') 



\x — x'|| 2 — (\r] — r]'\ — is)' 
r\r\' 



(30) 



(31) 



Note the distinction between the infinitesimal quantity e, used to define the pole prescription, 
and the parameter e = —H/H 2 . The propagator (I30p is the generalization of the Chernikov- 
Tagirov propagator for de Sitter space to space-times with constant, but arbitrary e |9J. The 
constant e propagator was already found for D = 4 in [TO] . 

Expression (1301) is precisely the result that was obtained previously by solving the prop- 
agator equation (j3J with the ansatz of (HH')^^ 1 times a function of y(x;x') [31 0]. If 
we employ the transformation formulae for hypergeometric functions and then their series 
expansion we can write this as, 



iAqq (x, x ) 



'l-e) 2 HH' 



(47T 



r(?-i)[. 



D 



-1 



a * 1 1 „ D y 



D-l 



fz/)r 



D-l 



ra+u)va-u) 



2-^1 



D - 1 



'l-e) 2 HH' 



(4tt; 



4n#-i 



r(2- 



2 

D - 1 



2 ' 4 



' 2 ' 4 



2 - 



T(\+u)T(\- V ) 



£ 

n=0 



r(|+V+n)r(|-V+n) (y^n-%+2 

T(3-f +n) (n+ 



/yy 



r(^ + z/+n)r(^-z/+n ) m n 
r(f+n)n! 



(32) 



(33) 



The problem with (I3"3"|) is that the gamma functions on the last line diverge for certain values 
of e, irrespective of whether or not x /fJ- = x M and with the dimensional regularization still in 
effect. For the inflationary case < e < 1 divergences occur at, 

2N „/£>-! 



r(— v + nj = T(-N + n) for JV = 0,1,2,... (34) 



D-2+2N 

For the decelerating case of 1 < e divergences are found at, 

e = 2^^-^Y(^ + u + n)=Y{-N + n) for iV = 0,l,2,... (35) 

3 Origin of the Problem 

The purpose of this section is to explain why the infinite space propagator ( 1331) diverges 
for the special values of e given in fl34H35p . We begin by discussing infrared divergences 
of the mode sum fl26|) . The discrete values of e result from how these infrared divergences 
are handled by dimensional regularization. Of course an infrared divergence should not be 
subtracted like an ultraviolet divergence! The correct procedure is to instead identify and 
remove whatever unphysical feature led to the infrared divergence. We close by doing this 
for the scalar propagator. 

It has long been known that the infinite space mode sum (|26|) has infrared divergences 
for broad ranges of constant e [11]. They follow from the small argument expansions of the 
Bessel function, 

Jv{z) = y i ±L W - , (36) 



and from its relation to the Hankel function 



HP(z) = Uz) + iN u {z) = —^-\ e-™J v {z) - J- U {z) . (37) 

The leading small z behavior for the terms on the first line of ( 1261) is universal, 

£=i T ( Ax \ z D ~ 2 ( Ax \^ 

z 2 x%( — z)^^x( — ) . (38) 



D-l , I LAX \ Z I LAX \ — 

2 \j^ Z ) ' r(^ li ) X \2y/rfrf 

What we get for the two Hankel functions on the second line depends upon whether the 
universe is accelerating or decelerating. For the inflationary case of < e < 1 the index v is 



7 



positive and it is the J_„ terms in the Hankel functions which make the leading contributions 
for small z, 

For the decelerating case of 1 < e < (D — 1) the index v is negative and it is the J +v terms 
in the Hankel functions which make the leading contributions for small z, 

'Mt = x ^ . (40) 




sin 2 (z/7r)r 2 (z/+l) 7r 

From the preceding discussion we see that the asymptotic small z form of the integrand 
in ( 1261) is a constant times, 

z D-2-2\u\ (41) 

For this to produce an infrared divergence requires the exponent to be —1 or less. For 
inflationary case of < e < 1 we have v > and the condition for a divergence is always 

1X16 1 

0<e<1 ^ _ 1 > fl _ 2 _^iz! = _ 1 _!(£z!). (42) 

1— e 1— e 

For the decelerating case of 1 < e < (D — 1) we have v < and the condition is met as long 
as e < 2(D-1)/D, 

l< e <?(£zi) ^ _ 1 > g -2 + ^£^-l- 2 ' J - 1 >-^. (43) 
D 1-e e-1 V ; 

Although the mode sum fl26|) has infrared divergences for all < e < 2(D — 1)/D, the 
final result (1331) only diverges for the discrete values of e given by expressions ( 131H351) . This 
is because dimensional regularization (T2J |T3] automatically subtracts power law divergences 
and only registers logarithmic divergences. For most values of e the infrared divergence is a 
power law, and dimensional regularization — quite incorrectly — sets it to zero. It is only 
for the discrete values (13411351) that a logarithmic divergence occurs and causes expression 
fl3"3"|) to become ill-defined. To see this, note that the logarithmic divergence could derive 
from any of the order z 2N series corrections to the leading small z term. For the inflationary 
range of < e < 1 the condition for a logarithmic infrared divergence corresponds precisely 
to flM}, 

0<e<l =^ (D-2) + 2N - D ~ 1 ~ € = -1 , (44) 

1 — e 

2N 

for AT = 0,1,2,... (45) 



D-2+2N 
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In D = 4 dimensions the problem values are e = 0, ~, |, |, . . . For the decelerating case of 
1 < e < 2(D — 1)/D the condition for a logarithmic infrared divergence agrees with f l35l ). 

-1 , (46) 
0,1,2,... (47) 

For D = 4 this corresponds to e = |, |, |, . . . 

It is important to understand that the infinite space propagator has physical problems 
for every value of e in the infrared divergent range < e < 2(D — 1)/ D, whether or not e 
happens to take one of the critical values (I34H35P necessary for a logarithmic divergence. This 
is because ultraviolet and infrared divergences mean different things. Ultraviolet divergences 
indicate that loop corrections have made an infinite change between observed parameters 
and the corresponding parameters of the Lagrangian. They can be canceled by expressing 
the parameters of the Lagrangian in terms of observed quantities plus counterterms which 
subtract off the divergences |14j . The automatic subtraction of dimensional regularization 
is not an error for ultraviolet divergences; it merely saves one the trouble of defining and 
subtracting the appropriate counterterm to cancel a power law divergence. 

Infrared divergences do not mean anything about parameters in the Lagrangian. Instead, 
they signify that there is something wildly unphysical about the computation being done. 
One does not deal with an infrared divergence by subtracting a counterterm; the correct 
procedure is rather to compute physically well-defined quantities. The classic example is 
the Bloch-Nordsieck switch from infrared divergent, exclusive processes to infrared finite, 
inclusive processes in quantum electrodynamics [2]. It will be seen that employing the au- 
tomatic subtraction of dimensional regularization to remove a power law infrared divergence 
corresponds to adding an illegal counterterm to make an unphysical question return a finite 
answer, instead of reformulating the question in more physical terms. 

The unphysical thing about the infinite space propagator is that a local observer cannot 
prepare the initially super-horizon modes of the state in coherent Bunch-Davies vacuum. 
Two plausible fixes have been proposed: 

• One could work on infinite space as in (TT3"|) but assume that the super-horizon modes are 
less singular than Bunch-Davies vacuum [T5|. Because only the super-horizon modes 
change there would be no effect on the Hadamard behavior of the propagator. If one 
continues to regard the state as obeying a(k)\Q) = 0, this fix corresponds to changing 
the super-horizon mode functions u(t, k) from the Bunch-Davies choice f|T5|) . Of course 
their time dependence is determined by the scalar field equation but their initial values 
and those of their first time derivatives can be freely specified. For example, if these 



1 < e < 



D 



(D-2) + 2N 



D-l- 

1-e 



2 D ~ 1 + N for N 
D + 2N 
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initial values were chosen to be those of the Bunch-Davies mode functions for e = D/2 
(regardless of the actual value of e) then there would be no infrared divergence, either 
initially or at any later time [16] . 

• One could also work on a compact spatial manifold such as a torus T D ~ X for which 
there are no initially super- horizon modes [T7]. In this case the free field expansion 
becomes a sum rather than an integral but it is generally valid to make the integral 
approximation to this sum, with a nonzero lower limit. When this was done for the 
graviton propagator on de Sitter background (e = 0) there is no disturbance to powerful 
consistency checks such as the one loop Ward identity [T8] and the nature of allowed 
counterterms [02 [20] ■ The renormalization of scalar field theories is not even affected 
at two loop order [2l"l 1221 123] . 

4 Finite Space Mode Sum 

The purpose of this section is to implement the second of the two fixes described above: the 
one based upon a finite-sized spatial manifold [17]. We show how this changes the mode 
sum for the propagator. We also derive the corrections it makes to the integrated, position- 
space form. Explicit demonstrations are given that the correction terms cure the N = 
and N — 1 divergences in expressions (I341l35p . And certain special cases are checked against 
known results [2Tj |24"] . 

We work on T D ~ 1 , which supports the spatially flat FRW geometry (Tj[|). If the coordinate 
radius in each direction is 27r//co then the integral approximation for the free field expansion 
of the operator is the same as (fT31 except that the integral is cut off at ||A;|| = k , 

<p(t,2) = f ( 2^i 0{k-k Q ){u{t,k)^ 3 a$) + t t *(t,fe)e~ i ^V(fc)} . (48) 

Of course the same cutoff works its way into the mode sum for the propagator f[26l) . 
iA(x; x) = 
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Here and subsequently z = kQ^Jrjrf. 

We can obviously break the integral over z up into two parts, 



dz 



2(1 



dz 



2(1 



dz 



(50) 



This means that the result for (T4"§1) is what we already have (153]) from the work of [21 El H] , 
minus the finite range integral. It would be simple enough to expand the integrand of this 
second contribution and then integrate termwise, but we really only need the most infrared 
singular parts. For the inflationary case of < e < 1 the index v is positive and the most 
infrared singular parts of the integrand derive from the J- u contributions to the Hankel 
functions. For the decelerating case of 1 < e < 2(D — 1)/D the index v is negative and it 
is the J +u parts of the Hankel functions that are the most infrared singular. We shall work 
out the series of leading corrections in each case. 

Let us begin with the most infrared singular correction for the inflationary case of < 
e < 1. From the small z forms (1381) and fl39|) we see that the desired correction is, 



5iA (x; x') 



:i- 



l HH' 



2 2 T 2 ' 



(4tt 



7T2r 



D-l- 



2(1 



dz z 



D-2-2v 



2 HH' 



--1 



T{2v)T{v) 2(1 -e) / 1 



(47T 



2(1-0 



(51) 



To reach the last form (15T1) we have used the doubling formula for the Gamma function 

2 2x - 



T(2x) 



7T2 



r x r 



i 

X+ 2 



(52) 



Of course the infrared divergence at z = in 5iA was dimensionally regulated, the same 
way as in the infinite space result iA^. This is wrong for SiA , just as it was for iA^, but 
expression (150!) implies that the two errors must cancel and we will shortly see this explicitly. 

We shall use the notation 5iA^ to indicate the N-th order correction in the case of an 
inflationary universe (e < 1). For the correction in the decelerating case (1 < e < 2(D—1)/D) 
we use the notation 5iA N . Before deriving the higher corrections, let us see that the addition 
of (1511) eliminates the N = divergence (1341) . Both flBTl) and the n = term from the last 
line of fl33|) have a common factor that we may as well omit, 



2 HH' 



--1 



1 



(47T 



X 



r(£+i/) 



(53) 
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The remaining contributions are, 

2(i- e )r(2z/)i» / i \^Ef r(f-i)r(2-f ) r(^i+z/)r 



2(l-e)r(2i/)r(i/) 



D-l 



r(|-^) 



r(f; 



e(£>-2)r( 



D-l' 



1 



<=(P-2) 



rf^i)r(i-f) r(^i+z/) r( z? - 1 



r(2i/) 



2(l-e) 
e(D-2) 



(54) 



It is convenient to expand this expression in terms of the small parameter a = e(D— 2)/[2(l— e)] 
such that in the limit when e vanishes expression (|54p reduces to 



2(l-e)r(2i/)r(i/) 



lim 

e-o e (D-2)r( 



D-l^ 



1 



£(£-2) 
2(1-0 



r(Z) — 1) < ln(aa') — 7r cot 



TlD 



+ 



2 In 



r 



D-l' 



D-l 



r(i-f)r v 2 



r( 



D-l 



i» r(i-iz) r(2«/) 



(55) 



Ho 



+ 1p 



D-l 



D 



2(1-0 
e(D-2) 

+ ^(D-l)-7^, 



where V , (- 2 ) = {d/dz)ln(T(z)) indicates the digamma function, we used a = —l/H Q r] (valid 
for e = 0), , 0(1) = —7 and the reflection formula for the digamma function, 



1p(l — x) = lf)(x) + 7T cot(7rx) 



(56) 



Multiplying (|55p by the common factor (|53p . and then adding the rest of f[3"3"j) — which is 
not singular for e = — gives the following result, 



limiAfx; x') 



H, 



D-2 



r(f) /4\f-i r(f+i) M\f-2 r(D-i 



t^0 



(4tt; 



-1 \y 



+ 



D 



-2 



+ 



ln(aa') 



™<f ) +2 Kf ) + *(^) - *(f ) +^>- i; 



n=l 



r(D-l+n) /y 



rxf+i+n) 

//Tf#+n) V4V ~ (2-f +n) (n+l)!V4 



n-f+2 



7 



O(fc 2 



(57) 



Except for the order k$ corrections, and for some constant, finite factors on the second line, 
expression ()5Tj) agrees precisely with the result first obtained in [21J and used subsequently 
in many one and two loop computations [TH [201 [22j [231 125] on de Sitter background. 
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It is straightforward to work out the next contributions from the lower limit. We merely 
add up the three first order corrections from the Bessel and Hankel functions for the case of 
v positive, 



Ax 

Jd-3 , Z 



"r/rj 



Ax \^ 



2y/rjrf 



D-3 

Z 2 



r( 



1 



Ax 2 
rjrj' 4 



' D- 



sin(z/7r)r(l — v) 
sin(z/7r)r(l — v) 



v 4 , n/J 



+ 0{z A ) 
0(z 4 



(58) 
(59) 
(60) 



The resulting lower limit term is, 

£-1 



(1- 



2 HH' 



2r(2z/)r(z/) 



(47T) 



r(|+i/)r( 



{rf+rj 



Ax 2 



4(^ — 1)7777' 2(D-l)r]r]' 



dz z 



D-2v 



(l-e) 2 HH' 



2r(2i/)r(i/) 



(47r)f 



r(|+i/)r( 



X 



[rf+r] 



/2\ 



Ax 2 



4(j/— 1)t7?7 / 2(79-1)7777' 



_£>+l-2v 
Z 

D + l-2i/ 



(61) 



This should cancel the N = 1 diverge of (I34p at e = 2/D, which affects the n = and n = 1 
terms on the last line of (EH, 



{l-e) 2 HH' 



£-1 

2 



r(^+*)r(£±i-i/)r(-§; 



2 



;^i) 2 -z/ 2 4 



(62) 



(4vr)f r(i+i/)r(|-i/) 

The key to seeing that the infrared divergence of flBTj) cancels that in fl62|) is to express 
both in terms of the small parameter, 

D+l\ (De-2) 



a = v — 



As before, we extract the common factor of, 



(l-e) 2 M' 



--1 



2(l- e ) 



(l-e) 2 M' 



(63) 



•—1 



(4tt; 



r(|+i/) 



(47T 



r(f +i+«) 



(64) 
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The lower limit contribution from (I6TT) is this factor times, 



2r(2i/)I» 



D-l> 



/2\ 



Ax 2 



4(^-1)7777' ~ 2(D- 1)7777' 



r(£> + l + 2a)r( 



fa 



D+l' 



4ar( , 

r(D+i) /2-y 



1+ 



„£>+l-2v 
~ z 

D + l-2t/ 



2 , /2 
77+77 

f]1]' 



2a 



As 2 



7777' 



1 \a 



V 



(2-1/) 



2^(^ + 1) +^ 



.D + l 



+ In 



77 2 + 7/ /2 
7777' 



^7777' 

+ 0(a 



In contrast, the contribution from (|62|) is (1641) times, 



r(^±i + ^)r(^i-^)r(-f; 
r(|-^) 



' 2 



'%l)2_j,2 4 



r(D+i+a)r(i-a)r(- 



4ar(-f -a 

rp+i) f /2-7/ 



(l + f)(l + a) 



a 



+ (2-2/) 



-V»p+l)+V(l)-^(- 



+2 



+ 0(a) 
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Adding (JHSD to (155]) and taking a = (which implies e = 2/1?) gives, 



T(Z? + 1) 



(2-y) 



K*sb) + ^ +1 > + *( 



Z? + l\ / £>\ 

"7- ^(-j J 



2 {rf+i 2 



D—l\ rjrj' 



rp+i) 



(2-y) 



D-2 



) ln(aa') +2 In 



(D-2)H 
Dk 



£+1 
D 



/"kD\ 

ncot[—) 



+^(D+1) + ~ 7 - V(y + l) 



/Z? 



2 /T] 2 + T]' 2 



D — l\ T]r\' 



+ 2 



We can get the full propagator for e = 2/Z? by multiplying flTOj) by the common 
S3]), and then adding the rest of (133]) with the now finite iV = correction fl5Tl) . 



lim iA(x; x') 



[(1-1) 2 HH'}* 



(4tt; 



r (f + 1 ) f^-\ r (f + 2 ) r 4 



(i-f)C-f) v ^ 



i r(f+3) 

'2! (3-f)(2-f) Vy 



r(z?+i) 
4r(f +i) 



(2-f)(l-f) V y 

2QP-1) 2 / 1 + t/ 
^777/ D — l\rj' 7] 



+2 + 1 + (2- v ){(l-i) ln(oa') -jrcot(^) + A' D } 



n=2 



r(f + 2+n) 



(2-f +n)(l-f+n) (n+l)!V4 



T(Z?+n) 



n(n-l)T(f +n) V4 



+ O(£; 2 ) . 



Here the constant i^Tp is, 



if D = 2 In 



+ ^(L>+1) + V(^) -7-^(y+l) 
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As far as we know the literature contains no result against which we can check ( 1711) but its 
limit in D = 4 dimensions has been worked out, 



.. .. . . , HH I 4 

iim urn iA i: x = - 

D^A t ^2 y ' > 64vr 2 



K " -2(i + *)-"»+i« 

I y k^rfg \r] g J 



+3(2-2/) [- \n(y) + 2 ln(||) + ± ln(aa') - 2 7 ] + 0(k { 



(73) 



This agrees perfectly with equation (3.82) of 

We have seen that the lower limit term which corrects the N = problem in (1341) is given 
by floTl) . For the N — 1 problem the corresponding lower limit correction is flBTl) . To see the 
general pattern, first substitute the relation for H in terms of g, 



H = H n 



-(l-e)H g 



This reveals the N = correction (lo"Tj) to be constant, 



(74) 



5iA 



:i- e ) 2 # 2 



(47T 



r(2i/)r(i/) 2(l-e) 

f r(i+f)r(^=i) ^=2) 



:i-^o 2 



2(1- E ) 



>2 1 / V 2 

The same substitution reveals that the N = 1 correction (IBTj) is quadratic, 



(75) 



5iAi = ± 



T{2v)T{v) 



D_ 

2 



(4tt)< M 2 



X 



4(i.-i) 2p-i; 



2(1-6) 
(De-2) 



Kq 



e(J3-2) 
2(l-<0 



(76) 



Both corrections are homogeneous solutions of the propagator equation 



d^fy^g^dJiA^ = 0. 



(77) 



Note that each lower limit correction Si An must separately solve fITTj) because each goes like 
a distinct power of ko- The freedom to add such homogeneous terms is precisely what is not 
fixed by just solving the propagator equation rather than using the mode sum. 
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We could work out the N-th lower limit correction SiA^ from the mode sum but that 
would involve tedious multiplications of corrections from the Bessel function and the two 
Hankel functions. A simpler technique is to use the fact that the correction must have the 
form, 



5iA 



-i — -l 



N 



r(2z,)i» 



2(1-6) 



(47T) 



rn+u)T(2=±) (D-2 + 2N)e-2N 



x 



,2 1 / V 2 
6(D-2) 



:i- e ) 2 ^o 2 



7.2 

«, 



2(1-0 



C E E *uA**rr U rf (S - k - ii ■ (7? 

fc=0 £=0 

Then we determine the coefficients by three requirements: 

1. The coefficient a^o derives entirely from the z 2N correction of the Bessel function and, 



(79) 



by direct examination of (j49j) . we can see that it is, 



jV!4^r(^l+JV) ' 

2. Symmetry under 77 -< — > 77' implies, 

Oh = ak(N-k-e) 

3. The series must of course solve (1771) . 
The differential equation (1771) implies, 

2i/-l 







<9 2 H d E E a ^ Ax 



2fc 2« /2(JV-fc- 



fc=0 £=0 



JV JV-Jfe 



2&(2A;+.D-3)a jW A 



k=0 1=0 

N N-k 



-EE 4£ ( £ - ^K^ 2 y*~v 



2fc^2£-2^/2(A r -fc-£) 



fc=0 £=0 
JV-1 N-l-k 



E 2(A; + l)(2A; + D-l)a fe+1 ^Ax 



2* 2< /2(tf-l-fc-/) 



fe=0 £=0 

iV-1 N-l-k 



E E 4(£+l)(£+l-7/)a fcm Ax 2 VV 



2fc^2£ /2(iV-l-fc- 



fc=0 £=0 



50) 



51) 



(82) 



(83) 



17 



Hence the coefficients must obey, 

{k + l)(2k+D-l)a k+u = 2{£+l)(£+l-v)a ki+1 . 
The unique solution consistent with the other two of the three properties is, 



1\N 

X 



W)r 2 (i- 



k\£\ (N-k-e)\ r(k+^)r(£+i-v)r(N-k-e+i-v) 

For N = this gives the known result , 

N = a 00 = 1 • 

A less trivial check is that it also works for N — 1, 

1 , 1 



N 



o-oo — a 01 



and 



a 10 



2(D-1) 



U) 



15) 



(86) 



17) 



Let us turn now to the decelerating case ofl < e < 2{D — 1)/D for which the infinite 
space propagator (1331) diverges at the discrete values given in (|35|) . By paralleling what we 
did for the inflationary case one can show that the lower limit term which corrects the N = 
problem is, 



(l-e) 2 M' 



■—l 



T{-2v)T{-v) 



(1 



(47t)t 
2 HH' 



dz z 



D-2+2u 



(47t)"? 



T(-2u)T{-v) 2(6-1) / 1 \ 
T(\-v)T{^) 2(D-l)-De\k 2 r ] r ] 'J 



2(D-1)-De 
2( e -l) 



(l-e) 2 # 2 



T{-2v)T{-v) 2(e-l) 



(4tt 



r(|-z/)r(^i) 2(£>-l)- J De 



(l-e) 2 ^o 2 



eCD-2) 
2(1-0 



2;; 



59) 



One can easily check that {k\r\r}')' lv solves the homogeneous equation ([77]) . So the full 
series of these lower limit corrections should take the form, 



5iA 



N 



(1-efHf 



T{-2u)T{-v) 



2(6-1) 



(l-6) 2 # 2 



A: 2 



r(|-z/)r(^i) 2p-i+iV)- (D+2iv)e 

2(l-e) JV /V-fc 

2k2£/2(N~k-£) 



(An j < ■ , o 

e(P-2) 



(90) 



fc=0 £=0 
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We determine the coefficients by the same three requirements as the a^, although the 
solution will be different because the ansatz (I9"0"j) is. 

We need to commute the differential operator in (1771) through the prefactor of {k 2 7]7]') 2L ' 
in the ansatz (190]) . 



d 2 + 



2v-\ 



V 



do 



kfov' 



2v 



1 2 / 

k rjT] 



2// 



d 2 + 



-2v-\ 
V 



(91) 



That is a highly significant result because it means the equation the bare series obeys is the 
same as we already solved for the lower series but with the replacement u — > —v. So we can 
write down the answer immediately, 



In N 



1 



r(^i)r 2 (i+z,) 



AJ k\£\(N-k-£)\ T(k + ^)T(e+l + u)T(N-k-e+l + u) 



(92) 



It is worth explicitly checking that the lowest N corrections 5iA N cancel the e poles in 
331) from the upper series ( 1351) . From ( 1901) and (1921) we see that the N = correction is, 



6iA° 



[{1-efHH 



2(e-r 



2(D-1)-De 
2(e-l) 



(4tt) 



_D/2 



r(^)r(|-i/) 



2(D-1)-De 



k\xiVj 



(93) 



This should cancel the divergence at e = 2(D — 1)/ D in the n = term on the last line of 

(USD, 

[(i-e) 2 M']f- 1 r(i-f)r(^+z/)r(^-z/) 



(4tt)^/ 2 

The relevant small parameter is, 

2{D-l)-De 



ra+u)va-u) 



D-l 



a = 



2(e-r 



(94) 



(95) 
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Adding (1931) to (I94p and taking the limit that a vanishes gives, 



lim 



[{l-efHH 1 }^' 1 T{^ + a)T{D-l + 2a) 



r(f+«)r( 



D-l- 



Q 



r( 



D-l> 



T(D-l+a) r(l- 



r(^i+a) r(L>-i+2a) r(i-f • 



-a 



r(l-a) 



[(l-i) 2 HH')^r(D-l) 



(4tt) d / 2 



r(f 



2 In 



(-Dtl-O-!)"^'' 



\ 2 J 



V 2 



D-l 



(96) 



(97) 



By using e — 1 = 1 — 2/D the final result for the propagator can be expressed in a form that 
is identical with the de Sitter case (IBj 



lim iA(x;x') 

, 2(g~l) 



[(l-e^M']^ 1 



r(f) M> r ' 



(4tt; 



;r(f+i) ^4 



f-2 



+ 



rp-i) 

r(f) 



(1 — e) ln(aa') 



-it cot ^ — — — J — |— 2 111 



h l^ 



D-l 



^(£)+v»(d-i)- 



"7 



n=l 



r(D-l+n) /y 



r(f +l+n) 



nr(f+n) V47 (2-f+n) (n + 1)! V4 



O(A; 2 ) . 



(9? 



Of course one must remember that in de Sitter e = and H = H , so (!57|) and (198]) are only 
formally the same. 

We will content ourselves with working out one more propagator. From (1901) and (1921) we 
see that the N = 1 correction is, 



[{1-tfHH']^- 1 2r(-z/)r(-2z/) / -1 

(4tt) d / 2 r(i-i/)r(^i) V4^v 



2Ax 2 r^ 2 +?7 



/2 



,2 \ 2 

where we define the small parameter a as, 

2D-(D + 2)e D+l 



D-l 



l + v 



2a 



(99) 



a 



2(c-l) 



(100) 
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This should cancel the divergences from the n = and n = 1 terms on the last line of 

[(i- e ) 2 M']T-ir(^±i+z/)r(i2±i-z/)r(-f^ r - D 



(47r) / 2 r(i+z/)r(i-i/) 

Adding ( |99i) to (11 01 1) and taking a to zero gives, 

[(l-e) 2 M']T-i r(£>+l) j 2 ^// // 



4 (^o 11 ) 2 -^ 



(101) 



(4tt; 



4r(f + 1) } D-l \rf + 7]) +2+ D 



"(2-y) 



L>-2 



) ln(aa') + 2 In 



D-2\H 
LVD + 2/ fc J 

£ + 1 



/Dtt 

71 COt , 

V 2 7 



(102) 



By taking advantage of the fact that (1 — e) = — (D — 2)/(D + 2) we can express the full 
propagator in a form identical to the N — 1 result ( !7TI) from the lower series, 



lim iA(x; x') 

[{l-efHH'}^- 1 \ T(f + 1) Mxf- 



2.D 
" D+2 



i r(f+3) /4\f-3 r(D+i) 



i r(f+2) /4\f-2 

+ (2-f)(l-f) ly 



2! (3-f)(2-f) Vy 



4T(f + 1) 



2P-1) 



n=2 



+2 + — + (2-y){(l-e) ln(aa') - vrcot( 
T(f + 2+n) 



1 + 1 
D — l\r]' i] 



(2-f +n)(l-f +n) (n+l)!V4 



T(D+n) 



n(n- l)T(f +n) V4 



Here the constant Cp is 
C D = 2 In 



k -I 



+ V(i>+1) + ^ 



D+l 



+ O(k 2 ) 



D 



(103) 



(104) 
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5 The scalar stress-energy tensor 

In this section we shall calculate the expectation value of the scalar stress-energy tensor 
using the propagator obtained in the previous section. The stress-energy tensor for a scalar 
field <f, with the Lagrangian (jHJ), is given by 



T, 



5S 



(105) 



The expectation value with respect to the vacuum state |0) can be written as 



, (106) 



where d'^ = ^J^. The propagator iA is given in terms of the infinite space propagator fl30l 



and the corrections flTHj) and fl90l . 

oo oo 

iA(x; a/) = iA OQ (x; x) + 5zAat(x; x) + (JiA^a;; a;') . 

N=0 

5.1 The infinite space contribution 



(107) 



N=0 



We first consider the contribution to fl 106[) coming from iA^. From (1311) we find that at 
coincidence (y — > 0) the following two identities hold 







2/=0 
j/=0 



-2(l-e)^V- 



(108) 



Moreover, since in dimensional regularization all .D dependent powers of y can be automati- 
cally subtracted, we find using (6.131.2) in [5] that the contributions from the hypergeometric 
function appearing in iA^ relevant for this calculation are 



_ (D — 1 D — 1 

d (D — 1 D-l 
2F1 — + v, — — 



dy 



D 


1 


v, — 


' 2 




D 




v\ — 


1 


' 2 





rq-f)r(f) 
y =o r(| + i/)r(i - u) 



y=0 



2D 



D-i\i\ r(i-f)r(f) 



r(| + i/)r(i - v) 



(109) 
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Using these identities we immediately find that 



dpffjAco (x; x' 



-H 



x 



D 



t r i^ r (i-f) r 

J. C r~> 



V) 



(47t)t 



r(| + „)r(± - v) 



D ~ 1 \ 2 2,0,0 1 ( 2 fD-l\2 



5V 



110) 



Making use of Eq. f II 6 j) the one-loop contribution to the stress-energy from iA^ can be 
written as, 



(4tt) d / 2 



r(i + v)r(|-v) 



L> - 1 
x I v 



D - 1 



o 

«-° + ^( o J •'//»' 



.(111) 



5.2 The e < 1 correction 

Next we consider the i5A^ contribution (1781) . We define for convenience 
An — — 



_h.2N 
ft 



(4vr 



k 2 



D-l 



v 



(112) 



such that 



N N-k 



5iA N (x; x) = A N y j y j a M (Axy if // 



2k 21 /2(N-k-e) 



(113) 

Since at coincidence Aa; is zero, we find that the only nonzero contribution to (11061) arises 
from k — 1 when both derivatives hit Ax, and from k = when both derivatives hit rj and 
rf. Thus at coincidence we have 



' N-l 



d p d' a 5iA N (x; x') 



2(N-l) 



N 



(114) 



2(iV-l) 



£=0 



where 



0x0 



Vp<? = Vpa + $ p S', 
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(115) 



Using the expression for ( 1851) we find that we can resum the series over I to obtain 
dp& a 5i/S.^{x; x) 



{-l) N T\l-v)T{N-\-v) } 
N v^4" T(iV)r(iV - v)T(N - 2v) 1 



(116) 



where we made use of the identity, 



N\ T(N+l-2u) = T(2N+l-2u) 



From Eqs. (11061) . (11121) and (11161) we find the following contribution to the stress-energy 
tensor from the iV-th correction, 

fO\T ION H °\ l ~e\ D (4) N+£ ^- v T(l-N + u)T(-N + 2u) 

WW)* ~ 2(47r)^/ 2 N+^-uT(I-N + u)T(^)T(N) 1 J 



(-(D-l) + 2v + (D-2)Ny 2 5y u + ^(-(D-l)-2(D-3)i/+2(D-2)iv) 



where we have transformed the Gamma functions such that they are not singular at the 
iV-th pole (14T1) and we used 

(Q-2)e 

H 2 (l ~ eff ~ l (^^j 2 " £ ' £ = ^11 " ^(^^-V , (119) 
with 2 = ka\r]\ at coincidence. 
5.3 The e > 1 correction 

Finally we consider the correction due to 5iA N given by Eq. (1901) . We define for convenience 
such that 

AT JV-fc 

^(x; x') = B N {ki m 'y Yl E M^) 2 yy 2(Ar - fe ^. (121) 

fc=0 ^=0 
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The calculation is similar as in the e < 1 case, and it yields: 



N-l 



d p d'JiA N (x;x') 



= B N (kWr[ ^6 W (-2^ 2(jv - 1 

e=o 

N 



(122) 



Summing over £ gives 



'-i) n a» r 2 (i + i/)r(iv - § 



2 1 ^ 2(JV-1) 

v^F t(n)t(n + i/)r(jv + 21/) ' 



D - 1 



(123) 



N 



resulting in the following contribution to the stress-energy tensor, 



(Q|T^|0> 



N 



2(4vr) D / 2 A^+^ + zy r(|-iV-z/)r(^±i)r(iV+l) 
^D-l)(l-2u)-(D-2)N)aX^ + \((D-l)(l-2u)-2(D-2)N)g, u 



(124) 



5.4 Renormalization 



The total one-loop stress energy tensor is the sum of the three contributions fillip . (11181) 
and (11241 . 

oo oo 

(QIT^Q) = (filr^lfi),*, + ^(niT^l^jv + £(^7^)" . (125) 

V=0 ^=0 

Note that the ultraviolet divergence (which in dimensional regularization appears as a term 
multiplying 1/(D — 4)) is confined to (II lip . Indeed, when expanded around D = 4, Eq. (II 111) 

gives, 

(niWc f {2 ~ e)eHD 1 



8tt 2 

(2-e)eH 4 
16V 2 



D-A 
4 - e + 7 + In 
3 



47T 

(2 - e)eff 4 
128vr 2 



1 



1-e 
^ + 0(D-4). 



1-e 



(126) 
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It is known that this theory can be renormalized by the R 2 counterterm only. Indeed, taking 
a functional derivative with respect to g^ v of the counterterm action results in 



-9 ^ 



J d D x^g~aR 2 = a{AV^V u R - Ag^DR + g^R 2 - ARR^) . 



(127) 



Making use of the corresponding expressions for R and R^ u in FLRW spaces (see e.g. Ref. [3]) 
this evaluates to 



S 



^Sg 



d D x^g~R 2 = 144(2 -e)e# 4 



-# 4 (48e(l - 4e + e 2 )a 2 5j5° - 12(3 - 22e + 22e 2 - Ae 3 )g^ (D-A) 



+ o((D-Af 

Note that the H D term in Eq. (1 1 2 6 [) can be expanded as, 



(128) 



1 + 



D-4 (H 2 

In — r 

2 V/x 2 



0{{D-A)< 



(129) 



where fi is an arbitrary renormalization scale. From Eqs. (11261) and (11281) we see that the 
divergence in (11261) is canceled by 



a 



D-4 



1152tt 2 (D-4) 



(130) 



where controls the undetermined finite part of a. The renormalized stress-energy tensor 
can be now easily obtained 



H 4 



1 



(2 + 16e - 16e 2 + 3e i )a 2 5°X - -(12 + 62e - 215e 2 + 146e 3 - 2Ae 4 )g i 



48vr 2 
3e(2 - e) ( 7 + V 



1 - e 



1 - e 



In 



1 - e) 2 H 2 
Ann* 



+ J2((n\T, u \n) N + (n\T, u \n) N ) 



N=0 



(131) 
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where the terms in the last line denote the D — > 4 limit of Eqs. flllSp and (I124p . 

The sums over N can be performed and the final result can than be recast in terms of a 
hypergeometric function 2F3. This is a useful procedure for studying the ultraviolet behavior 
of the corrections (when the cutoff \zq\ —>■ 00). Since here we are primarily interested in the 
infrared sector \zq\ <C 1, we shall not perform the summation over N. 



5.5 Resolving the divergencies of the digamma functions 

Even though the ultraviolet divergences have been removed by dimensional renormalization, 
the renormalized stress-energy (11311) still seems to diverge at the poles of the (di)gamma 
functions (l3"4T - [3"5l) (see also Eqs. fl4oTT4Tl) ). We shall now show that these divergences are 
only apparent however, and that they are canceled by the correction terms in (11311) given 
by (11181) and (11241) . precisely as there were designed to do. 

Let us first consider the case e < 1. Expanding (134)) around the iV-th pole we have 



-N + 5 



N-5 



v 



3 - + N 
2 



(132) 



1-e N + l-6' 

where N is a positive integer and 5 is an infinitesimal quantity. In this case the first digamma 
function in (11311) diverges, and its contribution to the stress-energy tensor is 

i£fwjH« + ( f -j)^] +0 ^- < 133 > 

where we used 

*P(-N + 6) = ~ + 0(8 ). (134) 


To check that our construction works next we rewrite the iV-th term (11181) by using fl 132[) 
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(135) 



When expanded in powers of 5 this gives, 



(Q|r^|Q) 



D^4 



N 



H 4 N(N+2) 
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^,Q_2V>(JV + 3)+V(3) 
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a ^S°5° + — 3 

N+l M 4(iV+ 1) 
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iV+l_ 
+ C?(5) .(136) 



1 One might think that the second digamma function in ()131|) has a simple pole also at e = 2, but that is 
canceled by the (2 — e) prefactor. 
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Note that close to the pole N/(N+1) = e + 0(5) and (JV-3)/[4(JV+l)] = e-(3/A) + 0(8), such 

are (to leading order in 8) identical. Moreover, 



that the tensor structures in 



and 



by comparing Eq. (I136P with (11331) we see that the 0(1/5) terms cancel, as required. The 
resulting leading order contribution to the one-loop stress-energy tensor is finite and depends 
logarithmically on the scale factor as oc iJ 4 a 2 [ln(a)-foonst.], where the constant term contains 
a logarithm of the cutoff k . 

In the case when e > 1 the singularities of the digamma functions in the stress-energy 
tensor f II 3 1 j) conform with the poles (|35|) . As above, expanding around the poles 



1 



-2-N + 5 



---N + S, 



N + 3-5 



1-e " 2 N+2-5 

we obtain the following contribution to (11311) from the digamma function: 

H A (iV+l)(JV + 3)l/ ,„n„n / 3^ 



(137) 



16tt 2 



(N + 2)' 



'138) 



On the other hand, the correction ( 11241) contributes the following to the scalar stress- 
energy (11311) 
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if 4 (iV+l)(iV+3) 
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N + Q 
4(iV+2) 



- 5 + ln(^ 2 ) -i!>(l)- 2^(iV+4) + ^(3) + 

3(5 / o rfl ff) 1 

2{NV2)\ a5 ^ + 2 9 ^ 



9»v 



N+2_ 
+ 0(5). (139) 



Noting that here (iV + 3)/(iV+2) = e + 0(5) and (iV + 6)/[4(JV + 2)] = e-(3/4) + 0{5), we 
see that the pole contributions 0(1/ S) in fl 139j) and (11381) cancel, resulting again in a finite 
contribution to the scalar one-loop stress-energy of the form, oc if 4 a 2 [ln(a)+ const.]. 



6 Discussion 



We studied the propagator of the massless, minimally coupled scalar on a spatially flat, 
FLRW background of arbitrary dimension and constant e = —H/H 2 . A previous result 
( 13"3"]) . first derived [H H] by solving the propagator equation was seen to agree with the 
Bunch-Davies mode sum for infinite space. It diverges at discrete values of e (I34ti35j) . even 
away from coincidence and with the dimensional regularization in effect. The origin of these 
divergences derives from the way dimensional regularization treats the infrared divergences 
that the mode sum possesses for all e in the range < e < 2(D — 1)/D. For most values 
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of e in this range the infrared divergences are of the power law type that are set to zero by 
the automatic subtraction of dimensional regularization. The infrared divergences become 
logarithmic for the special values of e given by the two series fl34l and fl35l) . which is why 
the infinite space propagator fl33|) diverges for these values. 

Of course one should never subtract infrared divergences, even power law ones. Infrared 
divergences signify an unphysical feature of whatever question is being posed and the correct 
way to avoid them is by making appropriate changes in the question. The unphysical feature 
of the infinite space propagator is that no local observer could prepare the initial state in 
coherent Bunch-Davies vacuum over more than a Hubble volume. A more realistic situation 
is attained either by preparing the initially super-horizon modes in some less singular state 
or else by simply not having any initially super-horizon modes. 

We implemented the latter fix, making the integral approximation to the discrete mode 
sum but with a nonzero lower limit. The resulting propagator (T491 could be expressed as the 
old result fl33|) . minus a series of homogeneous solutions. Canceling the infrared divergences 
requires only the most singular contributions. For the inflationary case of < e < 1 our 
results for the corrections are f!78|) and fl85|) ; for the decelerating case of 1 < e < 2(D — 1)/D 
they are f[9"U|) and f[9"2l . We showed explicitly that the first two divergences of (I54"|) and 
the first two divergences of fl5o^) are canceled in this way, and we obtained the resulting 
propagators with dimensional regularization still in effect. We also showed that these results 
agree with special cases which have been reported in the literature [2~Tl 124"] . Moreover as an 
example we calculated the one-loop expectation value of the stress-energy tensor and showed 
that the divergences as discussed in [I] are canceled and replaced with a term that grows as 
the logarithm of the infrared cutoff. 

The case of e = 1 deserves special comment. From the convergence of the lower sin- 
gularities fl34l) and the upper ones fl35l) it might seem that the infinite space mode sum is 
highly infrared divergent. However, this is just an illusion derived from the fact that our 
convention for the zero of conformal time shifts disco ntinuously at e = 1. If one takes the 
limit e — > 1 while holding fixed the co-moving time and the initial Hubble parameter then 
the mode functions have a perfectly well-behaved form, 
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The propagator also has a smooth limit, 




i 



(141) 
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where we define the symbol Y as 



Y(x; x) = H 2 Ax 2 - 



In 



- is) . (142) 



An important distinction exists between the inflationary case of < e < 1 and the 
case of deceleration. The reason for this is that the parameter used to characterize the 
infrared cutoff, zo, approaches zero in an inflationary space, but it grows without bound in 
a decelerating space-time. In an inflationary space-time the physical wavelength associated 
with the infrared cutoff, a(t)/ko, grows faster than the Hubble radius. Thus a super-Hubble 
cutoff stays super-Hubble at all times. In contrast, for a decelerating universe the Hubble 
radius grows faster than a(t)/k and therefore an initial super-Hubble cutoff will eventually 
enter the Hubble radius. Hence the effect of the cutoff becomes more and more profound 
as time evolves. This behavior can be seen nicely from the expectation value of the stress- 
energy tensor (11311) . The corrections from 5iA^ and 5iA N have the following respective time 
dependences: 

H\zl) N ~^ oc a ~* +N ^ 

H\zl) N+ ^ oc fl -*PH-"<i-«>] . (143) 

This behavior should be compared to the tree level contribution to the equation of motion, 
which follows from the Friedmann equations, and scales as H 2 oc a~ 2e and H oc a~ 2e . Thus 
we see that for an inflationary space, e < 1, all corrections N > 1 decay faster with time 
than the tree-level contributions, while the N = correction scales with time equally as if 2 
and H. For a decelerating space e > 1 however, there are infinitely many values of N for 
which both corrections grow indefinitely as the universe expands. @ 

Thus we conclude that for an inflationary universe, the initial conditions become less 
visible as time progresses. In this case it does not matter much which of the two fixes for 
the infrared problems of Section 3 is employed, or precisely how it is implemented, while the 
opposite is true for a decelerating universe. Although we have worked out the fix based on a 
finite spatial manifold, the absence of evidence for a finite size to the universe suggests that 
the more physically relevant fix for deceleration is the one based on a less singular ensemble 
of initially super-horizon modes [15] . Because the initial condition becomes progressively 
more visible for deceleration, the physically relevant choice is probably the one consistent 
with at least 60 e-foldings of primordial inflation. Note that modes which experience first 



2 In order to get the correct late time behavior in this case, one would have to sum the series in N and 
consider the asymptotic behavior of the resulting hypergeometric function 2-F3 in the limit when |zo| 3> 1- 
Since that limit corresponds to a sub-Hubble cutoff, it is not of great interest to us. 
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horizon crossing will be even more infrared singular than for Bunch-Davies vacuum. It would 
only be for currently unobservable super-horizon modes that one could use the less singular 
mode functions needed to regulate the infrared divergence. Hence one expects the infrared 
behavior of the physically relevant propagator for deceleration to be strongly influenced by 
the initial condition provided by primordial inflation. 
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